Theory of vortex structure in Josephson junctions with multiple tunneling channels: 
Vortex enlargement as a probe of ±5-wave superconductors 
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We theoretically study Josephson vortex structures in Josephson junctions which have multiple 
tunneling channels caused by multiple superconducting gaps. Deriving "coupled sine- Gordon equa- 
tions" from the free energy taking account of the multiple tunneling channels, we examine two typical 
situations, a heterotic junction composed of multigap superconductor, insulator, and single-gap su- 
perconductor, and a grain-boundary junction formed by two identical multigap superconductors. 
Then, we reveal in both situations that the magnetic field distribution of the Josephson vortex for 
=bs-wave superconductivity is more enlarged than that for s-wave without sign change between the 
order parameters. Its mechanism is ascribed to a cancellation of the multiple Josephson currents. 
We display such an anomalous Josephson vortex and suggest how to evaluate the enlargement. 

PACS numbers: 74.50.+r,74.20.Rp,85.25.Dq 



The discovery of iron-based high-Tc superconductor 
has triggered a numerous number of studies on its su- 
perconducting mechanism and properties. It is now 
well-known through various experiments—"— that multiple 
bands contribute to the superconductivity and multiple 
superconducting full (s-wave) gaps open below the transi- 
tion temperature. On the other hand, several theoretical 
works ^^^^ have proposed that a sign change occurs be- 
tween the 5- wave gaps when a strong repulsion works be- 
tween the quasiparticles on the disconnected Fermi sur- 
faces. The symmetry with such a sign change has been 
called ±5-wave— , and its peculiar features have been 
intensively exploredi^"— . 

In cuprate high-Tc superconductors, the experimental 
quest for the superconducting gap symmetry has a long 
history in which an epoch-making work was the detec- 
tion of a half quantized vortex in corner or tri-crystalline 
junctions The discovery was so conclusive that such 
a measurement has been regarded as the most reliable 
way to confirm unconventional pairing symmetry since 
then. Is such a type of phase sensitive measurement also 
available for identifying ±5- wave symmetry in iron-based 
superconductors? The answer is not so simple—, because 
it seems to be rather difficult for 5-wave case to detect 
the sign change in spatially twisted geometries. 

In this paper, we propose an alternative way based 
on the observation of a Josephson vortex to identiy ±5- 
wave symmetry. The size of the Josephson vortex unex- 
pectedly enlarges for ±5-wave compared to the size esti- 
mated without the sign change. Such an enlargement is 
widely observable in various junction configurations, e.g., 
a heterotic junction composed of an iron-based super- 
conductor, an insulator, and a single-gap superconductor 
(SIS) a grain-boundary junction formed by two iron- 
based superconductor grains—, an intrinsic Josephson 
junction only for highly anisotropic compounds-, and so 



on. The detection will be possible if one uses the scanning 
superconducting quantum interference device In this 
paper, we derive the "coupled sine- Gordon equations" for 
the Josephson junctions with multiple tunneling paths 
stemming from the multigap character. The equations 
predict an anomalous structure for the Josephson vortex 
in the ±5-wave case, in which the sign of the Josephson 
critical current density depends on the tunneling channel. 

The theory of Josephson junctions with multiple tun- 
neling channels is in great demand for examining and 
understanding weak link properties of multi-gap su- 
perconductors. A theoretical development was done 
by Brinkman et al^^ and Agterberg et a/.i^, as for 
MgB2 and NbSe2. The modification in the conventional 
Ambegaokar-Baratoff relation was shown in these lit- 
eratures. In addition, proximity effects were studied in 
a heterotic structure composed of a normal metal and 
a multuigap superconductor^^. The observation of col- 
lective modes in two-gap superconductors via Josephson 
junctions was also proposed We note that a peculiar 
effect of the sign change between the superconducting 
gaps on the Josephson current was suggested by Agter- 
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FIG. 1: Schematic for Josephson junctions with multiple tun- 
neling channels, (a) A heterotic junction between single- and 
two-gap superconductors, (b) A grain-boundary junction be- 
tween two-gap superconductors. 
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berg et al in a context irrelevant to iron-pnictide su- 
perconductors. After the discovery of iron-pnictide su- 
perconductors, the importance of examining the effects 
of such a sign change grows significantly. A large amount 
of studies about Josephson junctions or tunneling spec- 
troscopy have been reported, e.g., the Andreev bound 
states— —, the dc- Josephson effect and the Riedel 
anomaly—. On the other hand, a theoretical research 
about the magnetic properties of Josephson junctions 
with multiple tunneling channels has been never so far 
studied, except for our previous work^^. Thus, we de- 
velop theory of Josephson vortex in such a system on 
the basis of a microscopic approach for Josephson junc- 
tions— i^i. 

First, we examine the heterotic SIS junction The 
situation is shown in Fig.[Ha), where the electrode 2 (1) 
is a two- (single-) gap superconductor with width s (5'), 
and the superconducting phases are expressed as (/p^^^ and 
(^(2) ((^^). The free energy density on the zx plane— 
is given by 
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FIG. 2: (Color online) The single-vortex solution for the het- 
erotic junction. We set Lx — 2.5Aji, rj' — 10^, 771 = 10^, 
772 = 1.56 X 10^, 32/31 = 0.8, and |Jin|/ii = 5.0. (a) s-wave 
(Jin > 0). (b) ±s-wave (Jin < 0). (c) The magnetic field 
penetrating into the junction. The solid line is for ± s-wave, 
while the dash one is for s-wave. 
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Fj = -2^ — cosi9^'^ ^cosx, 



with respect to A^, and we have the coupled 

(2) sine- Gordon equations as 
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The label i represents the band index {i = 1,2). We de- 
fine and as, respectively, = {hc/e*)dx^^—Af and 
af = {hc/e*)dx(p^'^'^ — 7 where e* = 2e. The penetra- 
tion depth of the superconducting state in the electrode 
2 (1) is written as (A'). The third term in Eq. ([H 
describes the Josephson coupling energy. As shown in 
Eq. ([2]), the first term represents the contribution from 
two different tunneling channels, ji and ^2, and the last 
one is the internal Josephson coupling microscopically 
originating from an inter-band interaction The sign 
of Jin determines the relative phase difference, x between 
the two order parameters in the two-gap superconduct- 
ing electrode 2. When Jin < 0, the preferable value of 
X becomes tt, which corresponds to ±5-wave. The fi- 
nal term in Eq. ([1]) is the magnetic field energy, in which 
By = d-^{A^ - Af) - dxA\ 

Using Eq. (|3]) and the Euler-Lagrange equations with 
respect to A^, we have 



he 



(5) 



where rj' = A /s'd, rji = X'~ / sd^ and f]~ = Vi ~^ V2 - 
Equation ([5]) is called modified Josephson relation 
Combined Eq.(|5|) with the Euler-Lagrange equations 
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where AJ-^ = Airde^ji/hc^ and A^^^ = 47r(ie* | Jin|/^c^. 
The coefficients of the coupling terms have opposite signs. 
Then, we call the equations "± coupled sine-Gordon 
equations." We note that ± is not relevant to ±5- wave 
but common for Josephson junctions having multiple 
tunneling channels. 

Figure [2] displays a single Josephson vortex solution 
numerically obtained from Eqs. (|6|) and ([7]). The spatial 
scale is normalized by Aji, and the total length in the 
direction of x-axis is = 2.5Aji. The boundary con- 
dition is given by 0^^\-Lx/2) = 0, 0^^\Lx/2) = 27r, 
^(2)(_L,/2) = xo, and e^^\-Lx/2) = xo + 2^. The in- 
ternal phase difference x is automatically chosen so that 
the free energy becomes minimum. We initially choose 
— TT (0) for Xo when Jin < (> 0), and solve Eqs. ([6j) 
and ([7j) iteratively. For the junction parameters, we set 
T]' = 10^, r]i = 10^, r]2 = 1.56 x 10^, ^2/^1 = 0-8, and 
|^in|/ji = 5.0. The width of the current core of the 
Josephson vortex for ±5-wave as shown in Fig.[2fb) is 
much wider than that in Fig.[2fa) for 5-wave without 



3 



the sign change. Moreover, one finds an antisymmet- 
ric current pattern for ±s-wave^^. We can find that x 
is shghtly modulated around the vortex center, although 
it is almost fixed to be a specific constant (0 or tt). Us- 
ing Eq. ([5]), we evaluate the magnetic field distribution 
around the Josephson vortex, as is shown in Fig.[2fc). 
We find a significantly enlarged distribution for ±s-wave 
compared to 5- wave without the sign change. 

For further understanding of the above enlargement 
results, we turn back to Eqs. (|6]) and ([7]). When x is 
rigidly fixed as or tt, we have the following equation, 
which are asymptotically valid except for the Josephson 
vortex core (i.e., \x\ — > oo): 

^2^(1) ^ 7)i(x)sin^(i), dlO^^^ - 7)2 (x) sin ^^2)^ (g) 

where 7)1 (x) = (1 + + Vi)/><n + cosx (1 + ^0/'^j2 and 
mix) = (l + r/0A?i+cosx(l + r?' + r/2)/A^2 (x = 0, tt). 
Here, we emphasize that the characteristic spatial scale in 
Eq. dHJ strongly depends on the type of pairing symmetry. 
Since 7)i(7r) < 7)i(0), we claim that the solutions for ±5- 
wave are more widely spread than that for 5- wave without 
the sign change. This asymptotic analysis well explains 
Figs.[2](a) and[2](b). In addition, using Eq. (|5]), the mag- 
netic field distribution inside the junction is asymptoti- 
cally obtained as 



Lxd 
^0 



(9) 



where A = {l-\-r]'-\-f])~^Lx/27r. Thus, we clearly find that 
±5-wave leads to an enlargement of the magnetic field 
distribution due to fji{7r) < ^{0). The origin of such an 
enlargement is the cancellation between multiple tunnel- 
ing channels as shown in Fig.[2jb). In addition, we point 
out that the relative phase difference between the super- 
conducting gaps slightly fiuctuates around a fixed value 
when I Jin I ^ ji,j2— • The asymptotic forms are valid 
in this case. The above qualitative discussion, Eqs. ([8j) 
and (|9|) does not depend on precise values of the junc- 
tion parameters, as long as the condition is satisfied. On 
the other hand, a quantitative evaluation of the magnetic 
field distribution requires the detailed information of the 
junction parameters. We will discuss the quantitative 
way to identify the symmetry, i.e., 5- wave or ±s-wave at 
the end of this paper. 

Second, we examine the grain-boundary junction as 
schematically shown in Fig.lTJb). Both the electrodes 
are assumed to be identical (two-gap) superconductors. 
This type of junction is observed in a weak-link be- 
tween grains of a polycrystalline iron-based supercon- 
ductor— i^i^. Alternatively, the situation is theoretically 
equivalent to the intrinsic junctions stacked along the 
c-axis. The free energy density is basically similar to 
Eq. ([1]), but there are two differences. The first term in 
Eq. ([U) is substituted with — Y^'f=i{^/^^^i){^i,i)'^ ^ where 
afi = {he/ e'')dx^i^ ~ ^1- The Josephson coupling en- 



ergy term is replaced by 
2 
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where 6>2*i = (P2'' - - {e*d/hc)A2 1 and 



(p^^^^ . The first (second) term in Eq. (fTOj) is the intra- 
(inter-)band Josephson coupling energy between the two 
electrodes The inter-band Josephson coupling orig- 
inate microscopically from incoherent (momentum non- 
conserved) tunneling, which is the dominant process at 
rough boundaries. 

Repeating the same treatment as the previous case, we 
have the modified Josephson relation. 
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Another relative phase difference X2 is determined by the 
identity 



^2,1 ~ ^2,1 - X2 - Xl- 



(16) 



Equation ([T4j) can be regarded as the sine-Gordon equa- 
tion with respect to an interband phase difference 

Assuming the physical situation in which superconduc- 
tivity fully grows in each superconducting electrode and 
Xl and X2 are fixed as the same specific value inside the 
grain, we choose 6^2*1 (~^a^/2) = 0, ^2*1 (^^c/2) = 27r, 
and Xi(-^a^/2) = Xi(^x/2) = Xo as the boundary 
condition for the single Josephson vortex. We choose 
7r(0) as Xo when Jin < (> 0). From Eq. ([16l), we 
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FIG. 3: (Color online) The single- vortex solution for the 
Josephson junction between two-gap superconductors. The 
junction parameters are the same as in Fig.[2l except for 
312/31 — 311 1 31 — 0-6. (a) s-wave (Jin > 0). (b) ±s-wave 
(Jin < 0). (c) The magnetic field penetrating into the junc- 
tion. The solid line is for it s-wave, while the dash one is for 



have X2(-^a:/2) = X2{Lx/'^) = Xo- Figure [3] displays 
the single Josephson vortex solution. The ratios, ji, 
312/ ji = 32\l 3\ = 0.6, and the values of the other junc- 
tion parameters are the same as the previous. Figures 
[31(a) and [31(b) shows the shape of the single vortex solu- 
tion for 5- wave without the sign change and ±5- wave, 
respectively. The ±s-wave superconductivity leads to 
Xi = TT. Thus, Eq. ([16]) means that no phase difference 
between ^2^] and ^2^^ appears even though the electrodes 
are ±s-wave superconductors. We find that no anti- 
symmetric current pattern for ±5-wave appears. Figure 
[31(c) shows the spatial distribution of the magnetic field, 
which is evaluated via Eq. ([TT]) . We find the enlarging 
behavior for ±5-wave similar to the hetrotic SIS case. 
We also find the cancellation between the intra-band and 
the inter-band Josephson currents when Jin < (i.e., 
X\ = X2 = tt), from Eq. ([TQ|) . Thus, the magnetic field 
distribution of the Josephson vortex enlarges in the ±5- 
wave symmetry. 

Finally, we discuss how to detect the enlargement of 
the Josephson vortex size experimentally. Hereafter, we 
focus on the heterotic junction system. A similar discus- 
sion is possible for the grain-boundary junction. First, 
we determine a characteristic spatial length for s-wave 
without the sign change in the case of ji 32- The 
Ambegaokar-Baratoff relation then gives a simple for- 
mula. 
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FIG. 4: (Color online) The magnetic fleld penetrating into the 
heterotic junction, varying v — ^11,1 /rn,2. The spatial scale 
is Ac, which is evaluated via Eq. ([T7|) . We set Lx — 2.5Ac, 
ryi = 10^ r]2 = 2.0 x 10^ |Jin|/Jc = 10, and AVA^ = 0.82. 
(a) The approximate formula, Eq. Q for s-wave. The inset 
shows the comparison to the numerical solutions for s-wave 
(circle) when r — 1. (b) The comparison of the numerical 
solutions for ± s-wave to a s-wave vortex core size. 



erotic junction for the the 5- wave symmetry without sign 
change. The resistance of the junction in the normal 
state associated with the tunneling channel ji is written 



as Tn 



We note that R: 



-1 _ 



-1 

^n,l 



-1 
^n,2 



rn,iAVrn,2A^ We define A^ 



The quantity is the critical current density in the het- 



and 32/31 = 
as A^ = 2A^K{ki)/7r, 
where A^g = min{A%A^}, A[ = max{A%A^}, ki = 
[1 - (A|/Aj^)2]^/^ and K{ki) is the complete elliptic 
integral of the first kindi^^. The superconducting gap 
amplitude in the electrode 1 (2) is written as A^ (A*) 
as shown in Fig.lTJa). Here, an important point is that 
the quantity A* depends on only the superconducting 
gap amplitudes. The direct evaluation of each resistance, 
^n,i or rn,2, is not practical, but the combined one R^ is 
measurable in the normal state. Thus, one can evaluate 
from Rn and A^, both of which are supposed to be 
experim entally measu red, and define a spatial scale as 
= ^yhc^ /47rde*J^ . Normalizing Eqs. (|6j) and ([7]) via 
A^, we can find that the equations have a free parame- 
ter r = rn,i/rn,2- Next, we estimate the magnetic field 
distribution by employing Eq. ([9]). We also check the de- 
pendence of r on the magnetic field distribution. Figure 
m^a) shows distributions of the magnetic field obtained 
from Eq. (|9]) in various r for the 5- wave without the sign 
change. We then find that r dependence of the distri- 
bution is not at all significant. Hence, one can adopt 
the result for r = 1 as a theoretical prediction for no 
sign change. Figure IH^b) presents a comparison with the 
±5-wave case. We find that the field distributions in the 
±5- wave case are much wider than that in the case with- 
out sign change except for the cases, e.g., r = 1/9 or 
r = 9/1. When r = 1/9, j2/ji — 0.09, indicating that 
one of the multiple tunneling channels is inactive and the 
system is approximately described by a single-channel 
junction. This is not the case of the present iron-based 
superconductors. We emphasize that the magnetic field 
distribution of the Josephson vortex for the ±5-wave su- 
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perconductivity never obeys the prediction on the basis 
of the 5- wave without the sign change except for such ex- 
treme cases. In other words, the observed length of the 
magnetic field extent is much larger than the theoretical 
size for the s-wave symmetry. 

In conclusion, we studied the single Josephson vortex 
solutions in the heterotic SIS Josephson and the grain- 
boundary junctions, and revealed an anomalous enlarge- 
ment of the vortex core size for ±5- wave compared to 
the size estimated by the Ambegaokar-Baratoff relation 
for no sign change. All phenomena were explained on 
the basis of the cancellation between different tunnel- 
ing channels due to the ±s-wave superconductivity. As 
for the heterotic SIS Josephson junction, the cancellation 
appears between the two Josephson currents ji and j2- 



On the other hand, the cancellation between the intra- 
and inter-grain Josephson currents occurs in the case of 
the grain-boundary junction. Such a cancellation leads 
to an effective change in a characteristic spatial length 
(e.g., penetration depth). Consequently, the Josephson 
vortex widely provides a reliable way to detect the gap 
symmetry in iron-based superconductors. 
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